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Abstract 

A 'complete' framework for 77/7*7/7*7* interactions is presented. The 
emphasis is on providing a model for 77 physics at all photon virtualities, 
including the difficult transition region Q 2 ~ m 2 . 



77 physics is interesting in its own right, because of the challenges it offers to our 
understanding of complex QCD-related phenomena. In addition, 77 events may offer a 
significant background to other kinds of physics studies, such as SUSY searches. This note 
summarizes the model recently presented in [|I|]. It starts from the model for real photons 
in M, but further develops this model and extends it also to encompass the physics of 
virtual photons. The physics has been implemented in the Pythia generator ||, so that 
complete events can be studied under realistic conditions. 

A first step in a calculation is the flux of incoming bremsstrahlung photons. Here a 
machinery has been set up for the flux of transverse and longitudinal photons as a function 
of the photon virtualities Q\ and Q\ and their longitudinal momentum fractions y\ and 
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D2- Assuming isotropic azimuthal angles, W 2 can be calculated. The user can specify 
cuts on the range of these variables, in order to restrict the generation to interesting 
events. The flux is convoluted with 7*7* cross sections dependent on Q\, Q 2 , W 2 and the 
polarization states. Beamstrahlung involves only real photons, and is thereby considerably 
simpler. With the flux of such photons given by some external program, already the older 
machinery can cover this case. 

Photon interactions are complicated since the photon wave function contains so many 
components, each with its own interactions. To first approximation, it may be subdivided 
into a direct and a resolved part. (In higher orders, the two parts can mix, so one has to 
provide sensible physical separations between the two.) In the former the photon acts as a 
pointlike particle, while in the latter it fluctuates into hadronic states. These fluctuations 
are of 0(a em ), and so correspond to a small fraction of the photon wave function, but 
this is compensated by the bigger cross sections allowed in strong-interaction processes. 
For real photons therefore the resolved processes dominate the total cross section, while 
the pointlike ones take over for virtual photons. 

The fluctuations 7 — > qq (— > 7) can be characterized by the transverse momentum 
k± of the quarks, or alternatively by some mass scale m ~ 2A;_|_, with a spectrum of 
fluctuations oc dk 2 ± /k\. The low-fcj_ part cannot be calculated perturbatively, but is 
instead parameterized by experimentally determined couplings to the lowest-lying vector 
mesons, V = p°, u>°, <p° and J/ip, an ansatz called VMD for Vector Meson Dominance. 
Parton distributions are defined with a unit momentum sum rule within a fluctuation Q , 
giving rise to total hadronic cross sections, jet activity, multiple interactions and beam 
remnants as in hadronic interactions. In interactions with a hadron or another resolved 
photon, jet production occurs by typical parton-scattering processes such as qq' —> qq' or 

99 -> gg- 

States at larger k± are called GVMD or Generalized VMD, and their contributions 
to the parton distribution of the photon are called anomalous. Given a dividing line 
k ~ 0.5 GeV to VMD states, the anomalous parton distributions are perturbatively 
calculable. The total cross section of a state is not, however, since this involves aspects of 
soft physics and eikonalization of jet rates. Therefore an ansatz is chosen where the total 
cross section of a state scales like ky/k^, where the adjustable parameter ky « m p /2 for 
light quarks. The spectrum of GVMD states is taken to extend over a range ko < k± < ki, 
where k\ is identified with the p_i_ m in(s) cut-off of the perturbative jet spectrum in hadronic 
interactions, p±mi a (s) ~ 1.5 GeV at typical energies ||. Above that range, the states are 
assumed to be sufficiently weakly interacting that no eikonalization procedure is required, 
so that cross sections can be calculated perturbatively without any recourse to Pomeron 
phenomenology. There is some arbitrariness in that choice, and some simplifications are 
required in order to obtain a manageable description. 

A real direct photon in a 7p collision can interact with the parton content of the proton: 
19 ~ * 99 (QCD Compton) and jg — > qq (Boson Gluon Fusion). The p± in this collision is 
taken to exceed k\, in order to avoid double-counting with the interactions of the GVMD 
states. In 77, the equivalent situation is called single-resolved (or direct x resolved), where 
a direct photon interacts with the partonic component of the other, resolved photon. The 
77 direct (or direct x direct) process 77 — > qq has no correspondence in 7p. 

The space of 77 processes thus is three-dimensional, with axes given by the k±i, kj_2 
and p± scales. Here each kj_i is a measure of the virtuality of a fluctuation of a photon, and 
p± corresponds to the most virtual rung on the ladder between the two photons, possibly 
excepting the endpoint k ±i ones. So, to first approximation, the coordinates along the 
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k±i axes determine the characters of the interacting photons while p± determines the 
character of the interaction process. Double counting should be avoided by trying to 
impose a consistent classification. Thus, for instance, p± > kj_i with k±i < k and 
k Q < kj_2 < ki gives a hard interaction between a VMD and a GVMD photon, while 
k±i > p_l > kj_ 2 with k±i > ki and k± 2 < k is a single-resolved process (direct x VMD; 
with p± now in the parton distribution evolution). 

If the photon is virtual, it has a reduced probability to fluctuate into a vector meson 
state, and this state has a reduced interaction probability. This can be modelled by a 
traditional dipole factor {m 2 , / (m 2 / + Q 2 )) 2 for a photon of virtuality Q 2 , where my — > 2k± 
for a GVMD state. Putting it all together, the cross section of the GVMD sector then 
scales like 



*? dk 2 ± k 2 v ( Ak 2 



h\ ki ki \4k{ + Q 2 J 

For a virtual photon the DIS process 7*g — > q is also possible, but by gauge invariance 
its cross section must vanish in the limit Q 2 — ■> 0. At large Q 2 , the single-resolved 
processes can be considered as the 0(a 8 ) correction to the lowest-order DIS process, but 
the single-resolved ones survive for Q 2 — > 0. There is no unique prescription for a proper 
combination at all Q 2 , but we have attempted an approach that gives the proper limits 
and minimizes double-counting. For large Q 2 , the DIS 7*7 cross section is proportional 
to the structure function F 2 (x,Q 2 ) with the Bjorken x = Q 2 /(Q 2 + W 2 ). Since normal 
parton distribution parameterizations are frozen below some Qq scale and therefore do 
not obey the gauge invariance condition, an ad hoc factor (Q 2 /(Q 2 + m 2 )) 2 is introduced 
for the conversion from the parameterized F 2 (x, Q 2 ) to a (t^js' 



Q 2 Y Air 2 a em ^. 2 4n 2 a em Q' 2 
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Here m p is some non-perturbative hadronic mass parameter, for simplicity identified with 
the p mass. 

In order to avoid double-counting between DIS and single-resolved events, a require- 
ment p± > max(A;i, Q) is imposed on the latter. In the remaining DIS ones, denoted lowest 
order (LO) DIS, thus p±_ < Q. This would suggest a subdivision crL 7 DIS = cr^g — Ci-res> 
with Opjg given by eq. (||) and <JiJ res by the perturbative matrix elements. In the limit 
Q 2 — > 0, the DIS cross section is now constructed to vanish while the direct is not, so 
this would suggest 0l 7 dis < 0. However, here we expect the correct answer not to be a 
negative number but an exponentially suppressed one, by a Sudakov form factor. This 
modifies the cross section: 

( cr l-res\ 

^LODIS = a DIS — ^l-res > ^DIS eX P • (3) 



a 



7 7 
DIS 



Since we here are in a region where the DIS cross section is no longer the dominant one, 
this change of the total DIS cross section is not essential. 

The space of 7*7* processes is now five- dimensional: Qi, Q 2 , k±i, k± 2 and p±. As 
before, an effort is made to avoid double-counting, by having a unique classification of 
each region in the five- dimensional space. Some double-counting remain in the region of 
large x ~ Q 2 /(Q 2 + W 2 ), where an ad hoc suppression factor (1 — x) 3 is introduced for 
the resolved components. 



3 



In total, our ansatz for 7*7* interactions at all Q 2 contains 13 components: 9 when 
two VMD, GVMD or direct photons interact, as is already allowed for real photons, plus 
a further 4 where a 'DIS photon' from either side interacts with a VMD or GVMD one. 
With the label resolved used to denote VMD and GVMD, one can write 



7 7 



^DISxres, 



dir X res 



* * I ^resldir \ * * 

+ °rcsxDIS eX P I 7 * 7 * J + ^resxdir (^) 
\ ^resxDIS/ 



7*7* _ / W_ \ 7*-. 

^dirxdir + | Q 2 + Q2 + ^2 " ° 



Most of the 13 components in their turn have a complicated internal structure, as we have 
seen. 

An important note is that the Q 2 dependence of the DIS and direct photon interactions 
is implemented in the matrix element expressions, i.e. in processes such as 7*7* — > qq 
or 7*g — > qg the photon virtuality explicitly enters. This is different from VMD/GVMD, 
where dipole factors are used to reduce the total cross sections and the assumed flux of 
partons inside a virtual photon relative to those of a real one, but the matrix elements 
themselves contain no dependence on the virtuality either of the partons or of the photon 
itself. Typically results are obtained with the SaS ID PDF's for the virtual transverse 
photons m, since these are well matched to our framework, e.g. allowing a separation 
of the VMD and G VMD/anomalous components. Parton distributions of virtual longi- 
tudinal photons are by default given by some Q 2 -dependent factor times the transverse 
ones. The new set by Chyla |J| allows more precise modelling here, but first indications 
are that many studies will not be sensitive to the detailed shape. 

Some first studies with this new framework look promising || , but much further work 
remains to assess its usefulness. 
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